STUDY  NUMBER  C 


Circular  Error  Probability 
of  a 

Quantity  Affected  by  a  Bias 


-TrnTTfiirmT— 

■  ItfUl  aciBHrmc  ANB 

(wSfflCiL  iMFOWiATIOM 


Distribution  of  this  dooumont 
is  unlimited. 


June  1963 


\r~ 


U' 


USAF  Aeronautical  Chart  and  Information  Center 
Chart  Research  Division 
Geophysical  and  Space  Sciences  Branch 


CIRCULAR  ERROR  PROBABILITY  OP  A  QUANTITY 
AFFECTED  BY  A  BIAS 


Prepared  by 
Melvin  E.  Shultz 


Distribution  of  this  document 
is  unlimited. 


GEOPHYSICAL  STUDIES  SECTION 
GEOPHYSICAL  AND  SPACE  SCIENCES  BRANCH 
CHART  RESEARCH  DIVISION 
PRODUCTION  AND  DISTRIBUTION  PLANT 


NOTICES 


This  publication  was  prepared  for  Internal  use 
In  the  Geophysical  and  Space  Sciences  Branch 
and  does  not  constitute  an  approved  Aeronautical 
Chart  and  Information  Center  publication. 


i 


CONTENTS 


TITIE 
NOTICES 
LIST  OF  TABIES 
ABSTRACT 

1.  Purpose  and  Scope 

2.  Circular  Error 

3.  Elliptical  Error 
Rjj  -  Re  Differences 


5*  Conclusions 


liar  OP  TABIES 


TABU  WO. 
1 


TITLE 


PAGE 

21 


2 


Rc  "  x 


3 


k 


x 


22 

23 

2k 


111 


ABSTRACT 


A  procedure  for  determining  the  radius  of  the 
39*35#  50*  and  90$  probability  circles  for  a  biased 
distribution  Is  presented*  Both  circular  and  ellip¬ 
tical  normal  bivariate  distributions  are  considered. 
The  elliptical  distribution  Is  replaced  by  an  equiv¬ 
alent  circular  distribution  and  an  approximate  radius 
obtained  using  the  olroular  distribution  procedure. 
Tables  giving  an  Indication  of  the  error  resulting 
from  this  replacement  are  Included. 
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1.  Purpose  and  Scope 


In  error  analysis  It  Is  assumed  that  all  systematic  errors  have 
been  removed  from  the  observational  data.  This  is  seldom  true  for  It 
is  virtually  impossible  to  eliminate  all  systematic  errors.  In  the 
analysis  It  is  hoped  that  those  systematic  errors  of  consequence  have 
been  eliminated  leaving  only  numerous  small  systematic  errors  vhose 
combinations  cannot  be  distinguished  from  random  errors.  When  a  sys¬ 
tematic  error  of  consequence  has  not  been  removed,  the  value  obtained 
is  said  to  be  biased.  That  is,  the  value  deviates  from  Its  true  or 
accepted  value  by  some  Known  or  undetermined  amount.  Systematic  errors# 
their  origin,  form,  and  removal  are  discussed  in  ACIC  RP-2  and  AOIC  IR-96. 

This  report  is  concerned  with  the  effect  these  undetected#  biased 
quantities  have  on  the  probability  interpretation  applied  in  error 
analysis.  That  is,  it  Is  desired  to  determine  the  radlue  of  a  circle 
which  will  include  a  certain  portion  of  an  error  distribution  effected 
by  a  biased  quantity.  This  may  best  be  explained  in  the  following 
manner:  consider  a  missile,  with  circular  standard  error  ce  “  r#  aimed 
at  a  point  T,  (see  Fig.  l).  The  missile  is  biased  (or  the  point  mis- 
identifled)  by  the  amount  d;  therefore  the  distribution  of  Impacts  are 
centered  around  the  point  A  not  T.  It  is  desired  to  know  how  large 
the  Rc  circle  (centered  on  T)  must  be  to  Include  a  certain  proportion 
of  the  impacts  distributed  around  A.  Or,  given  limiting  values  for 
Rc  and  r,  how  large  a  value  of  d  can  be  allowed?  These  are  the  ques¬ 
tions  vlth  which  this  report  le  concerned.  Only  the  normal  bivariate 
error  distribution  will  be  considered. 


Pis.  1 

In  the  preceding  paragraph  the  error  distribution  is  a  special 
case,  that  is*  the  distribution  is  considered  circular.  To  have 
a  circular  distribution  the  standard  errors  in  the  x  and  y  directions 
must  be  equal  (<rx  *  Oy).  When  circular,  the  probability  function 
of  a  bivariate  normal  distribution  is  simplified  by  reducing  the  tvo 
parameters  crx  and  0y  to  one  parameter  oc*  Section  2  is  concerned 
vlth  the  determination  of  Rc  and  probability  values  under  the  cir¬ 
cular  condition. 

Unfortunately  it  is  unusual  in  bivariate  analysis  to  obtain  a 
circular  distribution.  The  general  case  is  an  elliptical  distri¬ 
bution  vhere  ox  y  Cy.  Since  this  complicates  the  solutions,  the 
elliptical  distribution  is  replaced  by  an  equivalent  circular 
distribution  in  Section  3  to  obtain  an  approximate  solution  by  the 
method  purposed  in  Section  2. 
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2.  Circular  Error.  If  no  bias  error  exists  (d  ■  0),  the  center  of 
the  error  circle  crc  coincides  with  the  center  of  the  probability 
circle  Rc.  Equation  (l)  expresses  the  probability  of  a  circle  of 
radius  Rc  with  a  circular  standard  error  oc. 

Pr  lx2  +  y2  S  R2]  -  [  exp  -  (x2  +  y2)  I  dxdy  (l) 

J  J  2n  oc  2oc  J 

■ 

Fig.  2  illustrates  equation  (l) : 

y 


Fig.  2 

When  the  center  of  the  ac  circle  does  not  coincide  with  the  Rc 
circle,  equation  (l)  takes  the  form: 

Pr  [x2  +  y2  £  R2]  -  J  J— exp  f-  K*  -  x)2  +  (y  -  y)2]ldxdy  (2a) 


3 


y 


pic.  3 

From  Fig.  3  it  can  be  seen  that,  if  the  coordinate  axes  are 
rotated  by  the  angle  a  so  that  d  lies  along  the  x-axls,  equation 
(2a)  can  be  expressed: 


Pr  [x2  +  y2  s  R2]  ■  f  f— i-g  exp  -  ~  [(x  -  d)2  +  y2]  dxdy 

J  J2it  oc  L  2 0% 


where :  d  ■  JT-  x)2  +  (y  -  y)2  for  all  coses. 


Fig.  4  illustrates  equation  (2b), 
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(2b) 


Fig.  k 

Convert log  equation  (2b)  Into  polar  coordinate* : 

x  ■  r  ooa  9 


*o  * 


y  ■  r  sin  6 


r-o  d-o 


R0  2* 


Pr  [Re]  -  exp 


['SW  j  <» 


r»»0  0«0 
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when  d  ■  0  equation  (  3»)  hM  t>»»  *?«§? 


Wben  A  y  0,  aquation  (3a)  cannot  too  Integrated  In  closed  form. 
Vltalls  (1956)  computed  his  table  of  Circular  Hemal  Probabilities 
toy  the  polynomial: 


Tabulated  valuea  of  probability  are  given  for  and  values*  The 

Re  * 

confutation  of  Pr  for  various  values  of  ~  and  ~  is  relatively  easy 

°c  Qe 


6 


with  a  machine  calculator  taut  very  time  consuming.  The  desired 
degree  of  accuracy  in  Pr  is  obtained  by  susmlng 

until  the  desired  correct  digit  le  unaffected. 

1 

0.60653 

Pr  -  0.6o653  [(1)  (l  -  O.60653)  +  0.5  (1  -  0.60653  •  1.5) 

+  0.125  (1  -  0.60653  •  1.625)  +  0.020633  (1  -  0.60653  •  1.645833) 

+  0.002604  {1  -  0.60653  •  1.646437)  +  0.000260  (l  -  0.60653  •  1.6)16697) 
+  0.000022  (l  -  0.60653  •  1.648719)] 

Pr  >•  0.60653  [0.44o4o6]  •  0.26712  or  26.712* 

The  table  was  utilised  to  develop  a  nomogram  and  graph,  figs.  5 
and  6,  as  an  uncomplicated  method  for  obtaining  values  of  Rg  and 
probability.  (The  probability  curves  of  39.35,  50,  and  90  per  oent 
were  graphed  because  they  are  the  probability  levels  of  circular 
standard,  CEP,  and  CMA8  errors.)  The  following  steps  should  be  fol¬ 
lowed  in  using  the  nomogram  and  probability  curves  t  enter  Pig.  5  with 

value;  place  a  straight  edge  at  the  origin  0  so  that  it  passes  over 
®c 

the  •$-  value;  obtain  numerical  value  from  the  nomogram  A.  With  the  A 

R 

value,  enter  Fig.  6  and  determine  an  — £  value  from  the  intersection 

°c 


n.' 

Example  1: 
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of  the  A  value  vlth  the  danirnrt  mmiyiMllty  curve.  The  final  Rq 
is  the  radius  of  the  probability  circle  resulting  from  a  normal 
circular  distribution  displaced  from  the  origin  by  distance  d. 


Example  2:  Given  a  missile  weapon  system  with  a  circular  standard 
error  (cc)  of  l?0m.  The  missile  la  aimed  at  a  point  200m  In  the 
x  direction,  and  4ocm  in  the  y  direction  from  the  target,  Vhat  is 
the  radius  of  the  $00  probability  circle  centered  on  the  target  T 


Assume  the  origin  of  coordinate  system;  (x,  y)  •  (0,  0)  as  the 
target. 


Solutions  _ _ 

d  *  yj{x  -  x)2  +  (y  -  y)2 


•c 


a.99 


enter  Pig.  5  with  y- 


^200)2  +  (400)2 
200,000  m  448a 


A  -  9.96 

enter  Pig.  6  with  A  -  9.96 
at  50%,  ?£  -  3.15 

c 


Rc  -  3.15  •  150 
-  472a 


The  result  (472m)  Is  the  radius  of  a  target  centered  olrole 
which  has  a  $0%  probability  of  containing  the  missile  Impact. 
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Third  order  annrrwliimt.lrm  for  39-35/  5-/  2nd  90^  prob¬ 

ability  levels  are: 


(§) 


$9.35 


-0-0498  (-g^)3  +  0.3085  (^)2  -  0.0076  (^)  +  1.010 


(?) 


50 


(S) 


90 


39.35 


±0.003 


-0.0535  ir)3  +  0.3952  (£)2  -0.0453  (£)  +  1.1B13 

0  GO 


'50 


±0.004 


-0.0550  (»-)3  +  0.3623  ($r)2  +  0.1674  (*:)  +  2.1272 


90 


-  ±0.011 


The  (§) 


39.35 


values  have  a  standard  deviation  of  ±0.003  with 
deviation  of  0.010  occurring  at  S  ■  0.  Values  of(^^  have  a  standard 
deviation  of  ±0.004  and  maximum  deviation  of  O.OOd  occurring  at 

-r-  •  3»00.  and  at  90jt  the  values  have  standard  deviation  of  ±0.011 

c  d 

vlth  a  saxlram  deviation  of  0.020  occurring  at  ■  0.45. 

ac 

Considering  example  2  vlth  ~  ■  2.99/  the  (3)^  value  vill  he 


(§>50  "  -°*°535  (  26.7306)  +  0.3952  (8.9400)  -0.0453  (  2.9900  +  1.1813) 


-  3.149  ±  o.oo4 

3*  Elliptical  Error.  As  stated  in  section  1/  a  bivariate  error 
analysis  generally  results  In  an  elliptical  error  distribution.  With 
Og  y  0y>  the  probability  integral  takes  the  form  of  equation  (4). 
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Zqu»lluu  (4)  is  not  integrapie  m  closed  form,  but  a  solution  is 

obtained  by  a  method  developed  by  Rosenthall  end  Hodden  (lyol)\ 

Evaluation  of  (4),  for  all  possible  values  of  (x  -  x),  (y  -  y), 

°mln  \ 

#•,  and  73-,  vould  result  In  a  table  of  monstrous  proportions, 
max  x 

therefore,  a  method  of  approximation  will  be  Introduced  so  that.  In 
many  Instances,  the  table  vill  not  be  required. 

If  d  is  zero,  x  and  y  are  zero  and  equation  4  Is  expressed: 

Pr  l*2  ♦  y2  *  •  J 

Equation  (5)  determines  the  radius  of  a  circle,  (x®  +  y^  s  R^),  vhich 
has  a  probability  density  equivalent  to  that  of  the  normal  elliptical 
distribution.  ACIC  TR-96  discusses  the  use  of  equation  (5)  for  re¬ 
placing  a  normal  standard  elliptical  distribution  with  a  more  convenient 
circular  distribution.  TO-96  suggests  that: 

a0  •  0.500  (at*  +  oy)  for  0.2  «  ^SiS  g  1.0 

J  WftX 

will  produce  a  standard  error  circle  (oe)  which  will  effectively  replace 
the  standard  error  ellipse.  By  replacing  the  standard  error  ellipse, 
the  evaluation  of  equation  (4)  Is  reduced  to  an  approximate  solution 
obtained  through  the  nomogram  and  probability  curves  or  equations  of 
section  2. 

1.  Ref.  1 


.4  U.£\ 

V*  'y  / 

« 


dxdy  (5) 
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Ex  ancle  3:  Given  a  mining  v<?e.pon  system  vith  a  crcasi'ange  cnw 

(ox)  Of  1.0  unit6  and  a  davnrange  error  (oy)  of  0.6  units.  The 
missile  Is  aimed  at  a  point  2.0  unite  In  the  x-direction  and  1.0 
unit  In  y-direction  from  the  true  position  of  the  target.  What 
la  the  radius  of  the  39.35$*  50$  and  90$  probability  circle  (Rc) 
centered  on  the  target?  Fig.  8. 

(x  -  x)  -  2.0;  (y  -  ?)  •  l.Oi  |  *  0.6 

oc  -  0.5  (1.0  +  0.6)  -  0.8 


enter 

enter 


d  -  7(2. 0)2  +  (1.0)2 
d  _  2.24 


2.24 


Fig.  5  vith  A 


0.8 

9-33 


2.80 


Fig.  6  for  39.35$  -  2.72 


Rc  -  2.1fl 


for  50$  -  2.98 


2.38 


l4 


4,  Rg  -  Rc  Differences.  An  error  in  either  probability  or 

P 

~  results  from  replacing  the  error  ellipse  by  the  corresponding 
circular  form.  This  error  is  the  difference  between  the  true  value 
of  the  elliptical  error  distribution  minus  the  approximation  value. 
A  br.’ef  discussion  of  these  differences  is  presented  in  order  to 
qualify  the  approximation  method  proposed  in  section  3*  The 
notation  in  this  discussion  will  be  as  follows:  The  Rc  value  Is 
the  radius  obtained  from  a  circular  distribution  by  equation  (3b). 
The  Rg  value  refers  to  the  radius  obtained  from  an  elliptical 
distribution.  Rg  will  be  used  in  general  cases  where  the  orien¬ 
tation  of  the  ellipse  is  not  in  question.  R?x  x  is  used  to  signify 
the  radius  obtained  by  equation  (4)  in  place  of  Rg  when  orientation 
is  considered.  The  first  subscript  denotes  the  direction  of  the 


larger  standard  error  value  and  the  second  the  larger  biased 
component,  or  R^  x  means  that  ox>Cy  and  (x  -  x)>(y  -  y).  The 
values  of  R0  are  those  obtained  from  VI tails'  table  and  the  Rg 
(R0  x,  etc)  values  are  obtained  from  the  tables  of  Rosenthal  and 
Rodden. 


For  the  39*35  and  50#  probability  levels  it  was  found  that 

if  the  ctppy  direction  and  the  direction  of  maximum  bias,  (x  -  £) 

or  (y  -  y),  coincide,  the  Rc  value  la  too  large.  This  problem  is 

shown  in  Fig.  9a  where  ffx»0y  "  *)>(y  "  ?)•  In  this 

particular  instance,  see  ex.  3,  for  50#  probability  RQ  ■  2.3d 

while  Rj  x  ■  2.35*  If  the  cn,v  direction  and  the  direction 
x, 
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of  maximum  bias  do  not  coincide  (Fig.  9b),  the  Rc  value  Is  too 
small,  xne  tt  in  this  situation  Is  2,41  units  and  the  R.  value 

v  ' 

is  the  seme  for  both  cases  2.33  units.  However,  for  some  of  the  90 $ 

probability  level  values  Ra^  y<  R0  <  Ra^  x. 

The  next  question:  with  the  same  error  ellipse  as  Fig.  9a  and 

9b,  what  are  the  R  and  probability  values  when  (x  -  x)  -  1,  and 

(y  -  y)  *>  2?  Since  d  Is  the  same  as  ex.  3/  the  RQ  value  for  50$ 

probability  level  is  again  2.3d  units,  Fig.  10.  Ry  rotating  Fig.  10a 

through  an  angle  of  90"  to  the  left  It  Is  a  reflection  of  Fig.  9a, 

thus  R_  „  *  R_  Similarly,  Fig.  10b  rotated  to  the  left  90° 

is  a  reflection  of  9t  and  R_  „  ■  R_  v. 

ffx,y  °y,x 

If  d  makes  an  angle  of  45*  with  the  coordinate  axis  (Fig*  11a 
and  b)  and  (x  -  x)  -  (y  -  y),  the  direction  of  ammv  is  Immaterial. 
Both  Fig.  11a  and  lib  have  equal  R  and  probability  values.  By 
rotating  11a  through  an  angle  of  90°  to  the  left  It  becomes  a 
reflection  of  lib.  Computing  Re  values  at  a  ■  45°  produces 
the  values  of  minimum. difference  (Rj.  -  Re),  As  d  progresses  in 
either  direction  away  from  o  ■  45s  the  difference  can  be  ex¬ 
pected  to  increase  to  maximums  at  0°  and  90° •  Tables  1  and  2 

give  values  of  (R0  v  -  Rc)  and  (Rc  -  R  )  respectively,  to 
x,y  °x. 

Indicate  a  possible  maximum  magnitude  for  the  differences.  The 

ft,,  j 

tables  cure  computed  for  the  - — "  ratios  of  0.3,  0.6,  0.4  and  0.2 

°max 

and  for  probability  levels  of  39*35#  50  and  90$.  Tables  3  and  4 
give  the  probability  difference  for  each  of  the  (Rg  -  Rc)  terms 
In  tables  1  and  2. 
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The  R  values  for  example. (3)  are: 


Rc 

39-35*  - 

2. Id 

Rc  x  ** 

°X,x 

2.09 

Rc 

50 

2.38 

R„  v  ** 
°x,x 

2.33 

Rc 

90 

3-33 

R-  „  ■ 
°x,x 

3.51 

From  table  2a,  b,  and  c  for  (x  -  X)  ■  2.0  and  (y  -  y)  ■  0,  the 
maximum  differences  to  be  expected  are  0.116,  0.66,  and  -0,189 
respectively.  The  values  In  the  example  all  fall  within  the  ex¬ 
pected  differences. 

5.  Conclusions.  The  simplified  computation  suggested  Is  feasible 
In  most  cases.  Where  the  error  In  the  R  value  created  by  the  trans¬ 
formation  from  elliptical  to  circular  forms  Is  great  enough,  such 
treatment  cannot  be  used.  Tables  1,  2,  3»  and  4  are  designed  to 
assist  in  determining  the  magnitude  of  such  errors.  Whether  to 
apply  the  less  complicated  method  or  not  Is  the  choice  of  the  user. 
The  Rosenthal  and  Rodden  method  must  be  used  where  the  error  Is 
considered  decisive. 
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For  39.35*  Probal 
(x  -  *),(y  -  y) 

>lUty 
(0.5,  0.0)  j 

(1.0,  0.0) 

(1.5,  0.0) 

(2.0,  0.0) 

Sin 

Sax 

0.8 

0.039 

0.062 

0.059 

0.6 

• 

0.089 

0.119 

0.1 18 

0.4 

0.149 

0.176 

0.164 

0.2 

m 

0.187 

0.209 

m 

For  50*  ProtebiU 
(x  -  x),  (y  -  y) 

ty 

(0.5,  0.0) 

(a) 

(1.0,  0.0) 

(1.5,  0.0) 

(2.0,  0.0) 

Sin 

Sax 

0.8 

0.033 

0.045 

0.036 

0.6 

m 

0.072 

0.076 

0.066 

0.4 

at 

0.107 

0.101 

0.082 

0.2 

- 

0.109 

0.103 

a* 

For  90*  PrcfcabiU 
(x  -  x),(y  -  y) 

^0.5,  0.0) 

(b) 

(1.0,  0.0) 

(1.5,  0.0) 

(2.0,  0.0) 

Sin 

Sax 

0.8 

-0.045 

-0.066 

-0.076 

0.6 

m 

-0.131 

-0.170 

-O.I69 

0.4 

m 

-0.255 

-0.28Q 

-0.342 

0.2 

m 

-0.388 

-0,423 

- 

(«) 

Table  2 
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rxuuttuiii ty  Differences 


For  39-35 1>  Probability 

(x  -»  x),(y  -  y)  |  (0.0,  0.5)  |(0.0,  1-0)  1(0-0,  1.5)  1  (0.0,  2.0) 


Probability  T11ff*r*ti£*5 


For  39*35^  Probat 
(x  -  x),(y  -  y) 

Rc 

Jllity 
(0.5,  0.0) 

(1.0,  0.0) 

(1.5,  0.0) 

(2.0,  0.0) 

°nln 

0BUUC 

0.8 

0.6 

0.4 

0.2 

For  50*  Probabil: 
lx  -  x),(y  -  y) 

m 

rty 

(0.5,  0.0) 

1.60# 

4.83 

7.73 

8.84 

(«) 

(1.0,  0.0) 

2.74* 

5.23 

7.39 

8.48 

(1.5,  0.0) 

2.4351 

2.71 

6.56 

m 

(2.0,  0.0) 

*Wn 

°max 

0.8 

m 

1.65* 

1.99# 

1.5U« 

0.6 

m 

3.65 

3.31 

2.75 

0.4 

m 

5.10 

4.21 

3.35 

0.2 

m 

4.86 

4.15 

m 

For  904  Probabili 
Jx  -  x),(y  -  y) 

ty 

(0.5,  0.0) 
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